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We investigate the scaling properties of phase transitions between survival and extinction (active- 
to-absorbing state phase transition, AAPT) in a model, that by itself belongs to the directed perco- 
lation (DP) universality class, interacting with a spatio-temporally fluctuating environment having 
its own non-trivial dynamics. We model the environment by (i) a randomly stirred fluid, governed 
by the Navier-Stokes (NS) equation, and (ii) a fluctuating surface, described either by the Kardar- 
Parisi-Zhang (KPZ) or the Edward- Wilkinson (EW) equations. We show, by using a one-loop 
perturbative field theoretic set up, that depending upon the spatial scaling of the variance of the 
external forces that drive the environment (i.e., the NS, KPZ or EW equations), the system may 
show weak or strong dynamic scaling at the critical point of active to absorbing state phase transi- 
tions. In the former case AAPT displays scaling belonging to the DP universality class, whereas in 
the latter case the universal behavior is different. 



I. INTRODUCTION 

The simple epidemic process with recovery or the Gribov process P, Q serves as an example of a prototypical 
nonequilibrium phase transition. It, also formally known as the Reggeon field theory is a stochastic multiparticle 

process that describes the essential features of local growth processes of populations in a uniform environment near 
their extinction threshold 0, 0] and belongs to the Directed Percolation (DP) universality class. For recent reviews 
see Refs. [H ||. A simple realization of the DP process is the predator-prey cellular automaton models. These 
models display nonequilibrium active to absorbing state phase transitions (AAPT) separating active from inactive 
or absorbing states due to competitions between spontaneous particle decay (death) and particle production (birth) 
processes. In these models for certain parameter values, the steady state density is zero, i.e., the species gets extinct. 
By definition, once reached an absorbing configuration, the system can not escape from these configurations [lfj |. 
Experimental realizations of DP universality has been reported in Ref. fill ] recently, where transitions between two 
topologically different turbulent states of nematic liquid crystals in their electrohydrodynamic convection regimes were 
observed. Ref. [ll| measured the relevant scaling exponents with high accuracy and found them to belong to the DP 
universality class. The critical behavior of the AAPT [8j depends on the conservation laws in the dynamics and the 
underlying symmetry. It has been conjectured in the form of the Directed Percolation Hypothesis [7], [l2| that in the 
absence of any special symmetry or conservation laws the AAPT belongs to the DP universality class as long as the 
system has a single absorbing state. In a more realistic situation, the spreading process in DP can also be long ranged. 
For example, consider infecting agents being advected by a local velocity field (e.g., parasites being carried in a wind 
flow in an ecological system) . To analyze such a situation, Ref. [l3| introduced a variation of the epidemic process with 
an infection probability distribution that decays with the distance r as a power law. Such long range DP models have 
analyzed in details numerically (l4j |. as well as analytically by using field theoretic methods [la, Il6j. In particular, 
Ref. [161 ] demonstrated four possible set of universal exponents for the long range DP problem, corresponding to four 
different pairs of (renormalization group) fixed point values for the two coupling constants in the model, reflecting 
several possible universal scaling behavior. Ref. [l6| enumerates the critical exponents at each of the fixed points and 
their stabilities. More recently, Ref. [17} discusses the effects of temporally (S-correlated multiplicative noises on the 
universal properties of the AAPT. They discussed the condition under which the usual DP universality class becomes 
unstable with respect to perturbations by the multiplicative noises. 

In the usual models belonging to the DP universality, the parameters defining the models are taken to be constants. 
Thus they represent a constant (non- fluctuating) environment. In this article we discuss the general question of how 
interactions with a spatially and temporally fluctuating environment affects the statistical properties of a percolating 
agent or the density of a near-extinct population <fi at the critical point of an AAPT belonging to the DP universality 
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class, by coupling the environmental fluctuations explicitly with the growth process. In order to address this question, 
we consider two different cases of fluctuating environments separately, namely, (i) a randomly stirred fluid, described 
by the Navier- Stokes (NS) equation, and (ii) a fluctuating surface, modeled by the Kardar-Parisi-Zhang (KPZ) or 
Edward- Wilkinson (EW) equations. We do not consider any feedback of <f> to the environment, i.e., the time evolution 
of the environment is autonomous. In all the cases, the coupling with the environment is nonlinear. Further we 
drive the environment (i.e., the NS, KPZ or EW equations) by Gaussian stochastic long-ranged forces with specified 
variances. Because of our choice for the dynamics of the environment, it itself displays universal spatial and dynamical 
scaling in the long wavelength limit, with the universality class depending upon the model used. In each of the cases 
we calculate the relevant scaling exponents at the critical point. The main result that we find is that the scaling 
behavior of the system near the extinction point of the AAPT depends generally upon the location of the system in 
the phase space spanned by dimensionality d and a parameter y that characterizes the spatial scaling of the noise 
correlations in the NS, KPZ or EW models, and depending upon the location in the phase space, the model may 
exhibit strong dynamical scaling, when the dynamics of the field <fi and the environment are characterized by the same 
dynamic exponent, or weak dynamical scaling, when the two have different dynamic exponents. Specifically in our 
model, weak dynamic scaling implies scaling properties of <f> belonging to the DP universality class, where as strong 
dynamic scaling corresponds to non-DP behavior. It may be noted that AAPT in contact with a randomly stirred 
velocity field has been considered in Ref. [l8j]. We compare our results and scheme of calculations with Ref. [l8| later. 

Motivations of our work are both theoretical and phenomenological. Our results are similar to those in Ref. [l6j . 
and provides for a mechanism to introduce long range flight in the otherwise local models of epidemic with recovery. 
Presence of both weak and strong dynamic scaling (in different regions of the phase space) makes it a good candidate 
to study general issues related to dynamical scaling in NESS. Apart from theoretical motivations, the models and 
results discussed here are useful in more biologically motivated context. Let us consider the specific example of a 
bacteria colony (or a biofilm of bacteria) undergoing simple cell division and death in its ordered (nematic or polar) 
phase. The nematic or polar order parameter being a broken symmetry mode obeys a dynamics that is scale invariant. 
This should be coupled to the dynamics of the bacteria density, which in addition to the growth decay terms, should 
have symmetry-allowed couplings with the order parameter field. Moreover, bacteria being living systems, there 
should be specific nonequilibrium stresses (or active stresses) arising from the continuous energy consumption of the 
living bacteria. Thus one obtains a coupled model of a density undergoing the extinction transition and a broken- 
symmetry order parameter field executing scale invariant dynamics. Our results here will help us understanding 
the interplay between the density field and the order parameter field in determining the universal scaling near the 
extinction threshold for more realistic naturally occurring systems. The rest of the paper is organized as follows: In 
Sec. Ill Al we discuss the basic DP model and the results from it in some details. Then in Sec. HTBl we discuss the 
case when the population density is coupled with a randomly stirred fluid described by the NS equation. Within a 
one-loop DRG approach we analyze the fixed points, calculate the associated scaling exponents determined by the 
stable fixed points and present a linear stability diagram. In Sec. Ill Cl we consider the environment to be a fluctuating 
surface, which is governed either by the KPZ or the EW equations. In all the cases we obtain the relevant scaling 
exponents by one-loop renormalized perturbative calculations. In Sec. Mil we summarize and discuss our results. 

II. EQUATIONS OF MOTION 

A. Directed Percolation model 

In order to set up the calculational background of our model, let us briefly revisit the problem of extinction of a 
single species in a uniform environment and its universal properties near the extinction threshold. As an example, let 
us consider population dynamics with population growth rate that depends linearly on the local density of the species 
and death rate that depends quadratically on the local density also undergoing a non-equilibrium active to absorbing 
state (i.e., species extinction) phase transition whose long distance large time properties are well-described by the DP 
universality class. In terms of a local particle number density <^(x, t) and taking diffusive modes into account, the 
Langevin equation that describes such a population dynamics is given by [see, e.g., Ref. Q] 

^ = DV 2 <i> + A 9 - \ d 4> 2 + v^C, (1) 

where D is the diffusion coefficient, \ g is the growth rate and Ad the decay rate of the density field (p. Stochastic 
function £(x, t) is a zero-mean, Gaussian distributed white noise with a variance given by 



(C(x,i)C(0,0)) =2D 2 S(x)5(t). 



(2) 
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The multiplicative nature of the effective noise in Eq. ([T]) ensures the in-principle existence of an absorbing state 
((f) = 0) in the system. Equation (TTJ) allows us to extract the characteristic length £ ~ y/D/\X g \ and diffusive time 
scale t c ~ £ 2 /D ~ D/\\ g \ on dimensional ground, both of which diverge upon approaching the critical point at A s = 0. 
Upon defining the critical exponents in the usual way 0] 

^(x,t->oo))~A|, (<j>(x,t))~t- a (\ g = 0), f~Aj w , *c~Q/X>~Ar^, (3) 
we identify the mean-field values 

[3 = l,a = 1,1/ = 1/2, and, 24, = 2. (4) 

Further, the anomalous dimension 77, which characterizes the spatial scaling of the two-point correlation function, 
is zero [9]. It is, however, well-known that fluctuations are important near the critical point and as a result mean- 
field values for the exponents are quantitatively inaccurate. In order to account for the fluctuation effects, dynamic 
renormalization group (DRG) calculations have been performed over an equivalent path integral description of the 
Langevin Eq. © @. A one- loop renormalized theory with systematic e-expansion, e = d c — d, where the upper critical 
dimension d c = 4 for this model, yields @, 

z = 2 - e/12,77 = e/12 and - = 2 + e/4. (5) 

The DP universality class, characterized by the scaling exponents (O above, is fairly robust, a feature formally known 
as the directed percolation (DP) hypothesis [12[. Only when one or more conditions of the DP hypothesis are violated, 
one finds new universal properties. For instance, the presence of long range interactions are known to modify the 
scaling behavior: Ref. [16| examines the competition between short and long ranged interactions, and identified four 
different possible phases. Our results assume importance in this backdrop: In our model, long-ranged interactions 
arise not because of any long range hopping, but due to coupling of the density <f> with an environment, modeled by 
a randomly stirred fluid or a growing and fluctuating surfaces, all of which in turn are driven by long range noises. 
It is expected that the presence of these long-range correlated background may affect the universal scaling behavior; 
our results below confirm this in general. 



B. Advection of simple epidemic process by a randomly stirred fluid 



Here we consider a randomly stirred fluid as a fluctuating environment coupled to the AAPT of a population density 
near its extinction threshold. 



1. Randomly stirred fluid model 



The dynamics of the fluid in the incompressible limit is described by a velocity field v which follows the Navier-Stokcs 
equation [l9[ 

dtVi + Ai(v ■ VH = D V V\ - — + fi, (6) 

P 

together with the incompressibility condition V • v = 0, which may be used to eliminate pressure in the usual way. 
Here p the pressure, p the density (a constant in the incompressible limit) and D v the kinematic viscosity. Parameter 
Ai is a non-linear coupling constant, which does not renormalize in the hydrodynamic limit under mode elimination 
due to the Galilean invariance of Eq. ([6]) (2(1 121| . Function f is the external force required to maintain a driven steady 
state. A promising starting point for theoretical/analytical studies on homogeneous and isotropic externally stirred 



zero-mean and Gaussian distributed with a variance 21 



fluid is the randomly forced Navier-Stokes model [20y2l|, where f is assumed to be a stochastic function which is 



2 V) 

(fi(^t)f j (-q,0))=S ij ^ = ^5(t) (7) 

in ci-dimensions, q is a Fourier wavevector, Di is a constant amplitude and y > 0. We are interested when d—A+y > 0. 
Since in this case the variance (O diverges in the hydrodynamic limit, the force f is said to be long-ranged and is 
infra red (IR) singular. This is a variant of the Model B (with the identification y — 2 — d) of Ref. [20( which 
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was subsequently used in Ref. [2l[ to calculate scaling behavior and various universal quantities associated with 
homogeneous and isotropic 3d fluid turbulence. The velocity field shows universal spatial and dynamical scaling 
independent of the microscopic viscosity and forcing amplitude [23]. These scaling properties are characterized by the 
roughness exponent Xv and dynamic exponent ze defined through the definition 

(<x,tH(o,o)>Hz| 2x */.(|z™, (8) 

where f v is a scaling function. Invariance under the Galilean transformation (see below) yields that the coupling 
constant does not receive any fluctuation correction in the hydrodynamic limit. A consequence of that is the exact 
relation between the exponents Xv and ze pU . I23I ]: Xv + %e — 1. Applications of one- loop DRG to systems with 
long range noises are more complicated and less controlled than its application in problems of equilibrium critical 
dynamics [2l|, |23|, |2J] . Despite the limitations, such calculations are successful in obtaining several useful results on 
dimensionless numbers and scaling exponents. Due to the infra red singular nature of the bare noise variance ([7]), 
the perturbation theory does not generate any fluctuation correction to the bare noise variance ([7]) which is more 
singular than it. Consequently, it is assumed to be unrenormalized. Thus in a renormalized perturbation theory, only 
kinematic viscosity D v undergoes nontrivial renormalization. An explicit one-loop DRG calculation yields 

z E = 2 - y/3. (9) 

Further, Xv = ~ 1 + 2//3. Thus for sufficiently high value of y, ze may even be less than unity, resulting into turbulent 
diffusion which is a very efficient way of mixing. In particular, the value y = d in (UJ is of particular physical interest, 
since it corresponds to the famous K41 energy spectrum for the velocity field: One finds for the one-dimensional 
energy spectrum E v (q) ~ k d ~ x (|i>j(q, t)\ 2 ) ~ g -5 / 3 in three dimensions. 



2. Extinction transition in contact with a randomly stirred fluid 



In the standard models for epidemic with recovery belonging to the DP universality class, the population density 
field 4> is allowed only to diffuse (apart from local reproduction and death). However in general, such processes (of 
reproduction and death) may take place in a fluid environment and the local population may get advected in addition 
to diffusing. Thus they may be called reaction- advection- diffusion systems. For instance, a reaction or birth/death of 
bacteria may take place in fluid, which in turn may be thermally fluctuating, or may be externally stirred. Here, we 
find out how the universal properties of standard AAPT are affected when the system is advected by an externally 
stirred fluid. 

The percolating field <p satisfies the same equation (TJJ), now supplemented by an advective non-linearity: 

^ + A 2 v = DV 2 <j> + \ g cf> - \ d cf> 2 + v^C, (10) 

Here, A2 is a coupling constant, which advectively couples <j) with v. Constants X g and Ad denote growth and decay 
rates. Gaussian noise £ has the same variance as ©. Redefining coefficients A g = Dt and A^ = -^y 2 - for calculational 
convenience, Eq. (TlOl) may be written as 

g + A 2 v • = D(V 2 + r)<P - ^p0 2 + v^C- (11) 

Thus the critical point is now defined by (renormalized) r = 0. 

As for the usual DP problem, the system exhibits a continuous phase transition from active to absorbing states as 
(renormalized or effective) r — > 0. The associated universal scaling exponents are formally defined as in © above. 
At the mean-field level the model Eqs. ([6]) and ([TTJ yield the same mean-field values for the scaling exponents as 
above in Sec lII Al Nonlinear couplings A 2 and Dg 2 , together with the expected large fluctuations near the critical 
point and multiplicative nature of the long ranged noise with variance (JTJ are expected to substantially alter the 
mean-field values of the exponents. Studies of these in-principle require full solution for the field $>(x, t). Equations 
© and (jlip being nonlinear, cannot be solved exactly. Hence perturbative means are necessary. We address this 
issue systematically via standard implementation of DRG procedure, based on a one-loop perturbative expansion in 
the coupling constants A 2 and Dg 2 about the linear theory. The resulting perturbative corrections to the correlation 
function may be equivalently viewed as arising from modifications (renormalization) of the parameters and fields in 
the model Eqs. © and (fTTj) . 

We being with the Janssen-De Dominicis dynamic generating functional [25j corresponding the Langevin equations 
© and (fTTj) together with the corresponding noise variances © and Q, which is given by 



Z NS = / DvDvD<j>D^eycp[S NS \, 



(12) 
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where (j> and v are auxiliary fields corresponding to the dynamical fields <j) an d v respectively which appear due to 
elimination of the noises from the generating functional -Zjvs- The action functional S^s is written as 

S NS (v,0,4>,4>) = Di y -0^ J dtv t v,P l ,(k)k 4 -y- d - J y dtvddtVi-i^Pijito^MtiviQt-tti + DukPvi} 

1 'dtfadt* - i A 2 A :i ^ Ui (q)^(k-q)+^(fc 2 -r)^+^^0(q)^(k-q)-^^0(q)0(k-q)}, (13) 



(2n) d J ^ ^ V ^ 2 ^rwry 2 

q q -i 

where Z?2 = A? s ) — — T^ 2 " i s the transverse momentum operator and Pjjj(k) = Py ;(k)fe; + P,/(k)fcj. Note 
that the last two non-linear terms in (| X3[) do not have the same coupling constant Dg/2, unlike the usual DP problem. 
This is consistent with the lack of invariance of Sns under the rapidity symmetry. The rapidity symmetry of the 
original DP problem (see, e.g., Q) is no longer admissible in the present case, since the Navier-Stokes Eq. §6$ being 
a viscous dissipative equation cannot be invariant under time inversion. 

Before we present our detailed DRG calculation, let us note the following: Since the dynamics of the randomly 
stirred fluid is independent of 4>, its dynamic exponent ze should be same as that in the absence of any percolating 
agent. We have ze = 2 — y/3 (see discussions above and the calculations below). Thus if there is a regime characterized 
by strong dynamic scaling, we should have z^ = ze = 2 — y/3. In this regime, the nonlinearity of the basic DP process 
(i.e., coupling constants g\ and 172) may or may not be relevant in an RG sense, corresponding to what we call LDP 
(long-range DP) and LR (long range) phases having different static scaling properties (the dynamic exponents are 
same in LDP and LR phases), characterized by the LDP and LR fixed points (FP), respectively. In contrast, when the 
coupling constant A2 that couples <j) and v is irrelevant in DRG sense, the dynamics of (f> is independent of v and hence 
<f> displays a dynamics that is identical to the usual DP problem with a dynamic exponent = zr>p = 2 — e/12. Can 
there be a phase displaying weak dynamic scaling with A 2 still being relevant? We expect not; because if A 2 is indeed 
relevant (in an RG sense), the the whole action (|13|) including the coupling term cj>v ■ V</> must be invariant under 
combined rescaling of space, time and fields characterized by a single set of exponents, i.e., a single dynamic exponent. 
Thus the assumption of weak dynamic scaling, i.e., the existence of two unequal dynamic exponents, rules out the 
DRG relevance of A2 . Finally one may in principle have a phase where all nonlinerities are irrelevant with the </>-field 
being characterized by the exponents of the linear theory (Gaussian phase) which turns out to be always linearly 
unstable. Thus, in short, we expect four different phases characterized by four FPs (LR, LDP, DP and Gaussian) and 
their associated set of exponents. Some of these phases may not be linearly stable. We shall confirm these physically 
inspired picture through detailed one-loop calculations below. In addition, we calculate all the critical exponents as 
well. 

Equations (j6|) and (|10|) are invariant under the Galilean invariance 

d d 

v^v + uo, 0^0, x^x + u f, t^t, — ->• — . (14) 

Invariance of the system under (|14p ensures that the coupling constants Xi and A 2 are equal and do not renormalize 
in the hydrodynamic limit: We henceforth set Ai = A2 = A below. The role of the (bare or unrenormalized) coupling 

constants in the ordinary perturbation theory in the present model is played by u — g\gi and w — —jyr"- In addition, 
there is a dimensionless number 8 = D/ D u , the Schmidt number which characterizes the ensuing NESS of the AAPT, 
and is a control parameter of the model. We set up a renormalized perturbative expansion in e = 4 — d and y up to the 
one-loop order. In order to ensure ultra-violet (UV) renormalization of the present model, we are required to render 
finite all the non-vanishing two- and three-point functions by introducing multiplicative renormalization constants. 
This procedure is well-documented in the literature, see, e.g., Ref. [26|]. Here, vertex functions of different orders are 
formally defined by appropriate functional derivatives of T[v, v, (p, <f>] with respect to various fields (dynamical and aux- 
iliary) , where T is the vertex generating functional and the Legendre transform of In Zns [HI ■ in the present model the 
following vertex functions show primitive divergence at the one-loop level: (i) / F. , (m)-t7tt, (m) r £, F r , . (iv) z £ F , . 



Their divergences may be absorbed by introducing renormalization Z- factors (see below). 

We employ the dimensional regularization scheme to compute the momentum integrals associated with the one-loop 
vertex function renormalization, and choose r = as our normalization point, where \i is an intrinsic momentum 
scale of the renormalized theory. The scaling behavior of the correlation or vertex functions in the hydrodynamic 
limit may be extracted by finding their dependence on /1 by using the renormalization group (RG) equation, which 
may in turn be obtained from the one-loop renormalization Z-factors. The Galilean invariance of the present model 
leads to exact Ward identities between certain two- and three-point vertex functions, which yields that the coupling 
constant A does not renormalize [27]. Further, since the bare noise variance (J7J is IR-singular, perturbation theories 
do not generate any correction to it which is more singular than it. Thus, the coefficient D\ also does not renormalize. 
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Fields v and v do not re normalize as well. The renormalized fields and parameters (denoted by a superscript R) are 
defined through the corresponding Z-factors as 

4> = Z(j} R , $ = Z0 R , t = Z t t r , D = Z D D R , g x = Z gi g R , g 2 = Z g2 g R , D v = Z Dv D R . (15) 

We perform explicit one-loop calculations in terms of coupling constants u — <?i<?2 and w = X 2 Di/D 3 / . From definition 
we have Z u = Z gi Z g2 and Z w = Z^ 3 . Now absorbing factors of 1/167T 2 into these coupling constants i.e, u —> u/l&ir 2 
and w — > w/l6ir 2 we can write down the Z-factors in terms of these scaled coupling constants. Further, since one of 
the Z-factors from the set Z, Z , Z gi , Z g2 , Zd, Z t , Zd„ is redundant, we use this freedom to set Z = Z without any 
loss of generality. We obtain 

Z Du = 1 - — (16) 

y 

Z = Z = l + ^, (17) 

_ « M -« \ 2 Di 1 2 f,-v 

Z ° 1 ~ 8 — ~ D„D{D U + D) {1 ~ d ] 16^ V ( } 

X 2 D 1 1, 2 u,-y uu- e 

D V D{D V + D) a 16tt z y 8 e 
, 3u/i" e A 2 L>i . l.it - * 

z " = ^YT^ w + D) ' 1 -^- (20) 

We hnd from the explicit one-loop results (|20|) that the Z-factors are linear in u, w and w = \ 2 Di/[D v D{D l/ + £))], 
(but not in 9 directly). Thus the perturbative expansions are in effect expansions in powers of u, w R and w. Further, 
Zq = Zd/Zd v is linear in w, but not in 9 itself. In contrast Z^ cannot be expressed as a linear function of w, along 
with u. In order to simplify the calculation, we make an approximation that 9 -C 1, such that w — 1/(D 2 D). This 
allows us to write Z^ linearly in terms of w. We show that, despite our simplifying assumption, we are able to obtain 
physically meaningful and interesting results which we present below. In the limit of small 9 we get Z-& = Z w Zd v / Zd- 
Using Z w = Zp , we can write Z^ = Z^ 2 Z^ 1 . Thus to derive Z^ we should find out what Zu is 9 -C 1. It turns out 
to be (after absorbing a factor of 1/167T 2 in the definition of w, setting d = 4 — e and keeping the lowest order term 
in e) 

u . 3w , . 

Zn = l-- M --- M -» (21) 

and hence 

z* = i + + — p-v + — fi-y, 22 

8e y Ay 

where a = §;jE§- The Wilson's flow functions can be defined as 

d d-d d d 

C<f = V-Q^ ln Z, Q = fi—\nZ, ( D = fi—\nZ D , ( T = n—\nZ T - 2, ( Diy = [i—ln Z D]j , (23) 



and the /3-functions as 



fa = n^-u R , fa = ^w R , fa = ^W R . (24) 



The renormalized coupling constants are written as u R = uZ u l fi e , w R — wZ^fi v , w R — wZ^ fi v , which for 
the /3-functions yields 

K - w R {-y + 3aw R ) , (3 U = u R (-e + ^u R + ^w R ^j , fa = w R (-y + 2w R + ^w R + ±u R ^j . (25) 

Fixed points of the model are given by the solutions of fa = fa = fa = which gives us various solutions depending on 
the different values of the the parameters u R , w R and w R . The only stable fixed point solution for wr is wr = yj (3a). 
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The different fixed point solutions for ur and wr are as follows: 

a ) u R — 0, w R — (the trivial Gaussian fixed point) (26) 
b ) u R = ^e,w R = (the DP fixed point) (27) 

c ) M fl = 0,w fl = ^y (LR, long range fixed point) (28) 

4 8 2 8 

d ) it fl = — e y, w R = e H y (LDP, long range DP fixed point). (29) 

5 15 15 15 

Since by construction u R and w R cannot be negative, from Eq. (|29[) it is obvious that the non-trivial fixed points may 
be present in the range Ay > e > |y. For e < |y the system has fixed points u R = and w R = |y which defines the 
LR fixed point. For e > Ay the system has fixed points u R — |e and w R = which defines the DP fixed point. We 
will see below that these ranges precisely coincide with the region of stability of those fixed points. 

The critical exponents are formally related to the Wilson flow functions [9j and hence to u R , w R and w R , and are 
given by 

u R 

V4> = V^ = -(* = -g~, (30) 
z E = 2-( Dv =2- aw R , (31) 

H = 2 -C D =2-^-^, (32) 

1 r „ 3u R 3w R , , 
v = -- = [*+ — + —]~\ (33) 

Thus the different exponents at different regions, defined by u R and w R , of the phase space spanned by e and y are 

• At the DP fixed point (^,0), the critical exponents are 

1* =, )r^'^ = 2_ ^' z E = 2- y/3, , i = 2 + |. (34) 

• At the LR fixed point (0, ^), the critical exponents turn out to be 

V4> = % = > z = Z E = 2 - | , ~ = 2 + |- ( 35 ) 

• At the LDP fixed point (4f — f§ , — f| + ff ), the critical exponents are 

m = % = ^-^, H = ZE = 1-\, I = 2 + i + |. (36) 

The interesting point to be noted here is that the dynamic exponent has the same value for both LR and LDP 
FPs, where has the corresponding static scaling exponents pick up different values. For example, at the LR FP the 
anomalous dimensions r]^ and 77^, which describe the spatial scaling of the correlation function and the propagator, 
are zero; hence spatial scaling of the correlation function and the propagator are given by the mean-field analysis (U), 
where as they pick non-trivial fluctuation corrections at the LDP FP. See Sec. IIIII for more discussions on this. 
To analyse the stability of the fixed points we must evaluate the matrix 

V du dw / 

and determine its eigenvalues A at u — u R and w — w R . The condition for infrared stable fixed points is that all the 
eigenvalues of the stability matrix M^s should be positive: 

• At the Gaussian fixed point (u R ,w R ) — (0,0) the eigenvalues are A = — e, — |. As both the eigenvalues are 
negative the fixed point is unstable. 

• At the fixed point (0, |y) the eigenvalues are A = — e + |y, The stability of this fixed point depends on the 

b 



condition e < %y. 
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For the DP fixed point (|e, 0) the eigenvalues are A = e, — + j^e. It is stable if e > Ay. 



For the fixed point u = 4e/5 — 8y/15 



,7,R — 



8y/15 — 2e/15, eigenvalues are 



A = 



10 5 



131e 2 32y 2 101ye\ 1/2 ' 



100 



75 



75 



/2 = Aj 



(say). Both eigenvalues are real. One of them A + is always positive within the window Ay > e > 2y/3 (this is 
the window in which this fixed point exists). Since detM^s = 2(3e — 2y)(y — e/4)/15 is negative outside this 
window, the second eigenvalue A_ must change sign when Ay — e or e = 2y/3. Thus this fixed point is stable 
for Ay > e > 2y/3. 

We see that the lines demarcating the regions of stability between the fixed points DP and LDP is given by w R = 0, 
which in turn is the same line where the values of the dynamic exponents z^ — 2 — e/12 and = 2 — y/3, corresponding 
respectively to the the DP and LDP fixed points, are equal, where as the line demarcating the regions of stability 
between the fixed points LDP and LR are determined by u R — 0. Further, as the system crosses over from the DP 
to LDP fixed point, z^ changes smoothly. The same is true for the anomalous dimension rj^ — rj^, which smoothly 
crosses over from its value e/12 at the DP fixed point to e/10 — y/15 at the LDP fixed point to at the LR fixed 
point. The correlation length exponent v shows similar behavior. In Fig. ([1]) below a phase diagram of the stable 
phases in the e — y plane is shown. Let us now compare with Ref. [l8T | where AAPT in contact with a randomly 







LR 






/ LDP 








DP 







FIG. 1: Schematic phase diagram depicting the stable phases of the model with a randomly stirred fluid as the environment. 
The continuous line which is given by Ay = e is the boundary between the LR and the LDP phases. The dashed line which is 
given by e = 2y/3 is boundary between the LDP and the DP phases. 

stirred fluid described the NS Eq. with a long-ranged force is considered within a onedoop renormalized perturbation 
theory like above. Our results for the scaling exponents and the phase diagram are same as that in Ref. 18]. The 
differences between Ref. [l8| and ours lie essentially in the details: Our choice for the coupling constants is slightly 
different from Ref. [llj]. Let us reconsider our choice for the effective coupling constants as used in the calculations 
above: The expressions of the Z-factors in ((20]) as well as Zg reveal that u, w and w = \ 2 Di/[D u D{D l , + D)] appear 
as the bare (dimensionless) expansion parameters in which the onedoop perturbative expansions are linear; but these 
are not linear in 9. Finally, our assumption of (bare) 9 <^ 1 allows us to write as a linear function of w as well. In 
contrast to us, Ref. [l8| worked with u, w and e (9 in our notation) as coupling constants in which the perturbative 
expansions are made. Our main operational motivation of expanding in terms of u, w and w is that u and w (or rather 
their renormalized counterparts ur and wr) directly describe the relative importance of the original DP nonlinearity 
vis-a-vis the advective nonlinearity. Hence, the plausibility of four possible phases described by {ur = = wr), 
(ur 7^ 0, wr = 0), (ur 7^ 0, wr 0) and (ur — 0, wr ^ 0) becomes immediately clear, a fact borne out by the 
detailed calculations described above. Nevertheless, there is no real contradiction between this work and Ref. [18[ as 
is evident by the same values for the scaling exponents and same phase diagram. 



C. Extinction transition in an environment of a fluctuating surface 



Consider next birth/death process of an agent taking place on a fluctuating surface: We thus now discuss the 
situation, where the dynamics of a population density field <f) near its extinction transition is assumed to be in contact 
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with a fluctuating/growing surface without any overhangs, represented by a scalar height field h(x,t), measured from 
an arbitrary substrate, which satisfies the nonlinear KPZ equation [3(| or the linear EW [28| equation of motion. 

1. Edward-Wilkinson surface growth equation 

The EW Eq.[28] is the simplest equation that describes a growing surface by a single valued height field h(x,t); 

dh 

— = D h V 2 h + vh. (38) 



We consider the noise tp is zero- mean, Gaussian distributed with a variance given by Eq. (1391) . We are concerned here 
with the effects of long ranges noises. The variance of the noise ip(lc, t) in the Fourier space is chosen to be 

(^(q, f)V(-q, 0)) = 2D iq 2 -y- d S(t), (39) 

where D\ is a constant setting the amplitude and the parameter y > 0. EW Eq. (|38[) is not invariant under the 
Galilean invariance. Eq. (|38[) . owing to its linearity, can be solved exactly. In particular, its dynamic exponent Zh = 2, 
regardless of the value of y. 

2. Kardar-Parisi- Zhang equation 
The KPZ equation [3(| is a nonlinear generalization of the EW Eq. above and is given by 

§ + y(V^) 2 -^V 2 / l + ^ (40) 

where h is the height field which gives the height of the growing surface from a reference plane, Ai a coupling constant, 
Dh is a diffusion constant, and tp is the external noise. Evidently, the EW Eq. (f3"5)) can be obtained from the KPZ 
Eq. (|40"]) by setting Ai = 0. Due to the Galilean (tilt) invariance (see below) of the KPZ equation the coupling constant 
Ai does not renormalize. Analogous to the NS Eq. ((5]) one defines dynamic exponent Zh and the roughness exponent 
Xh for characterization of the correlations of the fluctuations of h: One writes 

(Mx,i)Mo,o))~M 2Xh A(^ E A), (4i) 

where fh is a scaling functions. Nonrenormalization of Ai yields an exact relation Xh + ZE = 2. When if) is a zero-mean 
Gaussian distributed white noise, the KPZ equation describes a smooth to rough phase transition at d > 2 [29[. We 
will however be concerned here with the situation when the KPZ equation (|40[) is driven by a long range noise, same 
as (|39)) . Stochastic dynamics of the KPZ equation driven by a long ranged correlated noise has already been studied 
extensively in Ref. [31] by using DRG methods. Such applications suffer from several technical complications which 
are similar in nature to those for the NS Eq. © with a long-ranged noise. Nevertheless, one-loop DRG calculations 
have been successful is obtaining the scaling exponents. As for the randomly stirred fluid model, the only qua ntity 
that renormalizes here is the diffusion coefficient Dh- Explicit one-loop RG calculation yield ze = 2 — y/3 [3ll. |32|. 
an expression which is identical for a given y to that in the randomly stirred NS Eq. ^ . 

3. Extinction transition in contact with an Edward- Wilkinson fluctuating surface 

Let us now assume that the percolating agent <fr is coupled to a fluctuating surface described by the EW equation 
P8p. Thus an EW surface now forms the environment of the percolating process. On general symmetry ground the 
time evolution of the density field <f> may be written as 

^ + X 2 (Vh) ■ (V<f>) + \^V 2 h = DV 2 <f> + A g - \ d <p 2 + v^C, (42) 



where A2 and A3 are coupling constants. Equation (|42l) may be written as 

d(p 



nt + X(Vh) ■ {Vcj>) + A 3 V ■ (0V/i) = D^ 2 <f> + \ g cp - \ d 4> 2 + v^C, (43) 
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where A = A2 — A3. The A3-term in Eq. (|43[) . being a gradient term, is irrelevant in the long wavelength hydrodynamic 
limit (which is our region of interest here). Thus the effective equation for <f> that we are going to work with is 

^ + \{Vh) ■ (V4>) = DV 2 cf> + \ g <f> - \ d cf> 2 + V^C- (44) 

Just as in the previous case of a randomly stirred fluid environment, on grounds of general arguments we expect to 
find four different phases characterized by four different fixed points - Gaussian, DP, LDP and LR. The first two 
correspond to weak dynamic scaling where as the last two should display strong dynamic scaling. The EW Eq. 
being linear the corresponding dynamics has a dynamic exponent ze = 2. Thus, strong dynamic scaling for the 
density field <f> implies it will display dynamic scaling characterized by an exponent = 2. Although this corresponds 
to its value of the linearized theory [see the mean-field exponents given by Eq. , it does not really correspond to 
an AAPT characterized by the mean-field exponents, since all the other critical exponents (e.g., 77, v etc) differ from 
their mean-field values. Our calculations below confirms this picture. 

Using Eqs. and (J3SJ) we can write down the generating functional for the model given by Zew — 

J DhDhD(f>D(f)e~xjp[SEw}i where Sew is the action of the model. Function h and 4> are the conjugate auxiliary 
fields which appear due to the averaging over the noise distributions. Now to simplify the action and the calculations 
subsequently we rescale the fields as 4> i&4>i h ih and (f> — » j3<j) with a/3 = 1. The parameters are then rescaled 
as Dia 2 f5 = ^P-, A^a/3 2 = and A„ = tD. This modified action can be written as 



SEw(h,h,<t>,<p) = D 2 



jwfj ^^-J ^I h {d t h + D hq 2 h}-J ^ W -A^MqMk-q)q-(k-q) 
D 



+ Dk 2 4> -tD4>+- 5>!<Kq)0(k - q) - W(q)0(k - q)]}. (45) 
q 

Note that action (|45[) is not invariant under the rapidity transformation, similar to action (|13[) . As before, we are 
required to identify all the primitive divergent one- loop vertex functions. The vertex generating functional Tew 
is defined in the standard way as log Zew- The EW equation (|3"5)) being linear, there are no corrections to the 
vertex functions . and . where Tew is the Legendrc transformation of log Zew- However, there arc 

dhdh Shoh 

now non-zero one-loop fluctuation corrections to A. Thus, the vertex functions that are to be renormalized in order 
to render the action gSl) finite are (i) 5 r f f f , (ii) S ~ S <. E 7, , and (iii) ttttt 1 - The (bare) coupling constants for the 

o(po(p o<po(po(p O(po(po(p 

present problem are ur = gig 2 , w = -^3. In addition, dimensionless Schimdt number 9 = D^/D appears as a control 
parameter. Assuming renormalizability, we introduce renormalization Z-factors for each of the primitive divergent 
vertex functions: (i) ^-frrr , (ii) ttt^i (hi) tttW:, and (iv) frrfS ■ We use a minimal subtraction scheme, in which 
diverging parts of the associated one-loop fluctuation correction integrals are evaluated in inverse power series of 
e = 4 — d and y. 

The renormalized fields and the parameters are related to the corresponding bare quantities through the relation 

4> = Z R <t> R , j> = Z R 4> R , r = Z T r R , D = Z D D R , 9l = Z gi g R , g 2 = Z g2 g R and A = %A«. (46) 

We evaluate the Z-factors of the renormalized action in terms of the coupling constants u = g\g 2 and w = Di\ 2 /D 3 
and rescale these coupling constants with a factor of 1/167T 2 as we have done before in Section Hi Bl Writing Schmidt 
number 9 = ^ , the Z- factors are written as 



z = 


z = 




Z T = 


1 - 


3u 


Z D = 


1 - 


u 


z u = 


1 + 






1 + 


u 



U /i 

6 



(47) 



29{l + 9fy 
w(l - 9) 
26(1 + 9) 2 y 



»-y, (49) 



0(1 + 0)*/ 0(1+0)/ ' (5Uj 

(51) 



_^ ,,-v 



20(1 + 0)V 
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The Z- factor Zg can be easily calculated from Eq. (|49[) as there is no renormalization for D^: 
The renormalized coupling constants are defined as 

u R = uZ- 1 fi- £ , w R = wZ^-fx-y, , e R = ezg 1 , w R = wz-^r y . (53) 

The critical exponents are derived from these flow functions at the fixed points of the model which are evaluated 
easily from the beta functions 

Q Q O 

(3 U = [i—u R , &2 = fi—w R , pe = ^6 R . (54) 

OfJ, OfJ, OfJ 

At the RG fixed point then we have 

d / 5u R w R 3w R (l - 9 R ) 



8 e R {i + e R Y 26 R {i + e R ) 2 

a _ R ( _ 3 R w R (l - 6 R ) 2w R \ 

Pu ~ u { e+ 2 u 6 R {i + e R f e R (i + e R )J 



Fixed points of the model can be obtained from j3 u , and /3g by setting them all to zero. Note that the value of 
9 a at the fixed points cannot exceed 1. In fact the range is 1 < Or < 0. The different fixed point solutions are as 
follows: 

• Gaussian fixed point: u R = = w R , 9 R = 9 (bare value of the Schmidt number which is any finite positive 
number) . 

• DP fixed point:?^ = 2e/3, w R = 0, 9 R = 0. 

• Long range fixed point (LR): u R = 0, w R — 4y, 9 R = 1. 

• Long range DP fixed point (LDP): u R ^ 0, w ^ 0, 9 R < 1. We discuss this fixed point in details below. 
At the LDP fixed point values of u R , w R and 9 R satisfy the coupled equations 

3 n w R „ 1-9 R 

-> u ~ 2 nnn , qr\ ~fiW7~i , nm2 = e > ( 56 ) 



2 9 R (1 + 9 R ) 9 R (l + 9 R ) 2 

■> 9 R (1 + 9 R ) 2 9 R (1 + 9 R ) 2 
u R 1 — 9 R 

T-^ woro*r = =0 - (58) 



The solutions are given by 



. « _ 2 {by - 2e) 2 qR _ „ y 



e 



^ = 4(2y + e)/9, ri>* = 0* = 3f . (59) 

9 y - e 7y - e 



Fixed point values (|59| formally define the LDP fixed point. Positivity of 9r demands that solutions ([59]) arc 
meaningful only if either y > e or 7y < e, such that the ratio (y — e)/ (7y — e) is positive. However, when y < e 
and 7y > e, 9 is negative and hence solutions (I5T)1) are not physical. With the knowledge of the fixed point values of 
the coupling constants, one can now easily write down the corresponding scaling exponents: At the DP fixed point 
(^p 0, 0), the anomalous dimension T]^ = r/i, the dynamic exponent z and the correlation length exponent v are given 
by 

* = ** = ^>* = 2-^, ±=2+|. (60) 



12 



At the LR fixed point (0,4y, 1), the same critical exponents turn out to be 

ffo = Vj> = , z = 2 , i = 2. (61) 

Lastly, at the LDP fixed point we find 

'70 = ^ = ^( 2 ^ + e )> ( 62 ) 

i/" 1 = 2-7 T = 2+-(2y + e), (63) 

z = 2. (64) 

Thus, clearly both the LDP and LR fixed points correspond to strong dynamic scaling as the dynamic exponent of 
the density field z$ = 2 = ze, the dynamic exponent of the EW surface. Although this is identical to the mean-field 
value of z$, the scaling behaviors described by the LR and LDP fixed points do not correspond to mean- field scaling 
as can be easily seen from the values of the other exponents. 

Finally let us briefly consider linear stability analysis of the different phases, which can be performed in a standard 
way analogous to the previous case of randomly stirred fluid environment. We find, unsurprisingly, that the Gaussian 
FP is always unstable for any e > and y > 0. The DP FP is stable for e > 3y. In contrast the LR FP is always 
unstable, unlike the case when the environment is modeled by a randomly stirred fluid. We do not discuss the stability 
of the LDP FP due to the associated algebraic complications. However, the determinant of the stability matrix at 
the LDP FP vanish for y = e, suggesting borderline between stability and instability. Incidentally, the line y = e 
demarcates the region of existence of the LDP FP. Further, since the LDP FP does not exist for y < e < 7y and DP 
FP is unstable for e < 3y, there is no physically meaningful solution in the region y < e < 3y, a situation that does 
not arise when the environment is a randomly stirred fluid. Fig. (|2|) below shows the stable phases of the model with 
a fluctuating EW surface as the environment. 




FIG. 2: Schematic phase diagram depicting the stable phases of the model with an EW surface as the environment. The 
continuous line is given by y — e such that the LDP phase exists for y > e. The dashed straight line is given by y = e/3; the 
usual DP phase exists for e > 3y. There is no physically meaningful phase in the region between the lines y — e and y = e/3. 



4- AAPT in contact with a fluctuating KPZ surface 



We now briefly consider how a fluctuating KPZ surface may affect the universal properties of extinction transition 
of a population density tj> near its threshold. Due to associated algebraic complications our studies in this section 
are less extensive. Nevertheless we are still able to obtain physically interesting results consistent with the results 
obtained elsewhere in this paper. Field 4> follows Eq. (1441) . as in the case when the environment is modeled by a 
growing EW surface. However, unlike the case of the environment modeled by the EW Eq. (|38|) for the specific choice 
Ai = A2, Eqs. (|4U|) and (|4"2")l are invariant under the tilt transformation: 

h —> h + e ■ x, <t> <f>, x— s-x — Xiet. (65) 

As previously, Eq. (|42p reduces to Eq. (14"4"|) after discarding total derivative terms in the long wavelength limit. 
Function £ is a zero mean Gaussian white noise with a variance @. As before, the multiplicative nature of noise 
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ensures that the absorbing state <j> = is a solution of Eq. (l44l) . Similar to the case of a randomly stirred fluid 
environment, on general physical ground we expect four different DRG fixed points to exist, corresponding to four 
distinct phases: (i) Gaussian, (ii) DP, (iii) LDP, (iv) LR. The Gaussian FP is expected to be always unstable. The 
DP FP corresponds to weak dynamic scaling, where as phases corresponding to the LDP and LR FPs should display 
strong dynamic scaling. 

For performing DRG calculations, we use a path integral formulation in terms of the Janssen-De Dominicis gener- 
ating functional corresponding to the KPZ Eq. (|40p together with the Gaussian long range noise with a variance (|3U)) 
for correlation functions as before. Rescaling as for the case with a fluctuating EW surface as the environment, the 
action functional now reads [after dropping total spatial derivative terms, see the discussions preceding (|45|) ] 



SKPz{h,h,(j>,< 



-- 2D 2 J J dthhk 2 -y- d -J -0L J dth{d t h-^-^2 q -(k~ q )h(ci)h(k- q ) + D h k 2 h} 

- J 7^ J d4[d t 4> - \J2 q • (k - q)Mq)<Mk - q) + Dk 2 <p - Dt<P 

+ yEMq)^-q)-ft^(k-q)]]- (66) 

q 

Here again h and <fi are auxiliary (conjugate) fields. Like (|13j). action (|66[) is no longer invariant under the rapidity 
symmetry. Consequently, nonlinear coefficients g\ and g 2 in general are unequal. As for the usual DP problem, the 
system exhibits a continuous phase transition from active to absorbing states as (renormalized or effective) r — ¥ 0. The 
associated universal scaling exponents are formally defined as in above. At the mean-field level, the model Eqs. P0|) 
and (fl^tl) yield the same values for the scaling exponents as in Sec. Ill Al Just like the model in Sec. Ill B[ nonlinear 
couplings A and Dgi, Dg 2l together with the expected large fluctuations near the critical point and multiplicative 
nature of the noise with the long-ranged variance substantially alter the mean-field values of the exponents, a fact 
which we confirm below by our one-loop DRG calculation. Although in general g\ ^ g 2 , it is only the product .91,92 
that appears in the perturbative expansion. 

The structure of the perturbation theory and its renormalization is very similar to those in Sec. Ill Bl above. Due 
to the long ranged nature of the noise variance (|39j) there is no renormalization of D±. As in the previous case in 
Sec. Ill C 31 the role of the coupling constants in the ordinary perturbation theory in the present model is played by 

1 A 2 

u = g\g 2 , w = -Tjg-, w — -pg-.We again consider e > 0, y > 0, for which non-trivial critical exponents will ensue. The 

h 

vertex generating functional Tk pz is defined as the Legendre transformation of log Zkpz ■ Galilean invariance of the 
action functional (1551) ensures that the three-point vertex function S F' CI ' Z does not renormalize in the hydrodynamic 

' ' ShohSh " ^ 

limit. The vertex functions, which must be renormalized in order to render the present model renormalized, are (i) 
a2 f A 'f z , (ii) 52 ffA z , (iii) TTrrrr , (iv) Krrnf, and 5 4rrnT ■ Using standard methods as described in the previous 
Sections we in details perform one-loop multiplicative renormalization by introduction of the renormalization Z- 
factors, which render the theory UV finite. We use dimensional regularization together with minimal expansion to 
enumerate the renormalization Z-factors. From the renormalization Z-factors, one may then derive the RG flow 
equation in the usual way. Finally, the scaling behavior of the correlation or vertex functions may be extracted by 
finding their dependence on fj, by using the RG equation derived below. The renormallization Z-factors for the fields 
and the parameters are defined as 

4> = Z<f> R : , $= Z4> R , r - Z T r R , D = Z D D R , g x = Z gi g R , , g 2 = Z g2 g R , A = Z{\ R , and D h = Z D ,D R . (67) 

We evaluate the Z- factors of the renormalized action in terms of the coupling constants u = g\g 2 and w = D 2 /Df l 
and rescale these coupling constants with a factor of 1/167T 2 as we have done before in Section Til B 21 The different 
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Z-factors are 
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3w 9-1 
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1 + 36 1 
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where b = y/4d^. 

The renormalized coupling constants are written as 

u R = uZ- 1 ^ , w R = wZ- x pT v 



w R = w 



z~^-y, e R = 9z-\ 



(68) 
(69) 
(70) 
(71) 
(72) 
(73) 
(74) 

(75) 



The critical exponents are obtained from the Wilson's flow functions at the fixed points of the model which are 
evaluated easily from the zeros of the beta functions 



d d d d 

Pu = fJ-^-u R , /3 W = fi^—w R , ,8,3 = n—w R , (3 g = fi^—9 R . 
Ofi o/j, o/j, Ofj, 



(76) 



At the RG fixed point then we have 
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= w R 


Pu 


= u R 
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= w R 
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= 9 R 



-y- 



,R 



2mV 



3uV 



9 R (1 + 9 R ) 2 26 R {1 + 9 R ) 1 + 9 R 
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2w 
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e R {l + 9 R ) 9 R {1 + 9 R )9 R + 1 



1-9 R 



26 R (1 + 9 R ) 1 + 9 R 



bw J 



(77) 



It is clear that w = y/36 is the only stable fixed point solution for w . The other solution w — is unstable 
always and we ignore it from our discussions below. While we find the FPs and the corresponding scaling exponents 
below, we do not discuss their linear stability. Although the latter can in principle be done just as we do for the 
other models, it is algebraically much more complicated due to the structure of the fixed point equations. As before, 
we expect four different fixed points to exist, similar to the previous cases. We also expect weak and strong dynamic 
scaling in different situations. We find 

• u R — 0, w R — (Gaussian fixed point). The corresponding scaling exponents are those of the linearized system. 



,R _ 2 C -r _ 



e,w = DP fixed point. The exponents are given by those of the DP universality class. 



u R = 0; w R and 9 R may be solved from the coupled nonlinear equations 

w R 1 - 9 R 



29 R {1 + 9 R )1 + 9 R 



-bw R = m/3, 



2m; 



3uV 



1 



9 R (1 + 9 R ) 2 26 R {1 + 9 R ) 1 + 9 R 



2v w H v w H - 



V9 (1 + 9 R ) 2 



y- 



(78) 
(79) 
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Explicit solutions of the above equations are a difficult task, owing to their highly nonlinear nature. However, 
without their explicit solutions, one may already obtain the following information: (i) Since w R > 0, one has 
9 R > 1, thus D > Dh, (ii) Dynamic exponent = 2 — -yj, = 2 — 2 e R {i+e R ) i+e R = 2 — y/3 = ze and hence 
strong dynamic scaling, (iii) correlation length exponent v is given by v~ l =2 — 7* = 2 + y/3, and (iv) the 
anomalous dimension rf^ = n*^ = u R /8 = at the LR FP. 

• LDP fixed point: u R 7^ 0. Actual enumeration of the fixed point values of the coupling constants are very 
difficult due to the complicated nonlinear structures of the underlying equations, and will not be discussed here. 
Using, however, the approximation 6 3> 1, we obtain u R = 4(e + 2y/3)/7. Fixed point values 9 R and w R are 
to be obtained from the coupled nonlinear equations 

(80) 

6 R {l + 9 R ) 2 — — ^ e -R il + eR)2 -- ( 81 ) 

With w R = y/3b the above two equations may in principle be solved and solutions be obtained with the 
overall approximation of large 8 R . We do not solve these here explicitly. Nevertheless, we can already extract 
useful information without having to solve for the coupling constants explicitly. We find: (i) dynamic exponent 

z^ — 2 — 7^ = 2— ^ 9g«(i+e fl ) i+e R = ^ ~ Z//3 = ze and hence strong dynamic scaling, (ii) correlation 

length exponent v is given by i/" 1 = 2 - 7* = 2 + 3u R /8 - w R (l - 6 R )/[26 R (l + 6 R ) 2 ] =2 + u R /A + bw R = 
2 + (e + 2y/3)/7 + y/3, (iii) anomalous dimension r]^ = r]'~ = u R /8 = (e + 2y/3)/14, and (iv) by using positivity 

of 9 R and ui R , Ay > e from Eq. (|80|) for physically meaningful solution. 

Thus, with a KPZ surface as a fluctuating environment for an extinction transition that, otherwise (i.e., with a uniform 
environment) belongs to the DP universality class, the broad emerging picture is similar to the other two models of 
fluctuating environment considered here before. One generally finds both weak and strong dynamic scaling in different 
regions of the phase space spanned by e and y. The details, including the values of the scaling exponents, may of 
course depend upon the actual model of the fluctuating environment. Lastly, some technical comments regarding 
alternatives to the DRG procedure here is in order: As we commented before, the one-loop DRG procedure for the 
KPZ Eq. (14TJ1) with long ranged noise suffer well-known technical problems. As an alternative to it, self-consistent 
mode coupling method (SCMC) [33j and functional renormalization group (FRG) |34| have been used to extract 
large length-scale, long-time physics of the KPZ Eq. with long ranged noise. The SCMC expansion, as illustrated 
in Ref. [3j|, yields results which match with those in Ref. [31| at Id, where as for d > 1 the results of Ref. [33[ 
differ substantially from Ref. [31j and yields much more physically sensible results in the limit of short ranged noise. 
Similarly, Ref. [34| uses the well-known Cole-Hopf transformation and applies the FRG (up to two-loop) on the 
resulting partition function to obtain scaling exponents for the KPZ Eq. (|40[) with long-ranged noises. Their results 
clearly highlight the short comings of the one-loop DRG procedure. While the qualitative picture that emerges out 
of our one-loop DRG calculations here are expected to remain on general physical grounds, it will be interesting to 
calculate the details of the AAPT transition in contact with a KPZ surface by using the SCMC or FRG methods as 
illustrated in Refs. [13, The EW Eq. ([38} being linear such technical issues as for the KPZ Eq. do not arise. 
Nevertheless, investigation of the associated AAPT transition (in contact with an EW surface) by using SCMC or 
FRG methods would be useful. 



III. SUMMARY AND OUTLOOK 



This article is a study of how non-trivial fluctuating spatio-temporal dynamics of the environment may affect 
the usual directed percolation process with constant environment. We have separately considered cases when the 
environment is a (i) randomly stirred fluid described by the Navier-Stokes equation with a long-ranged force, (ii) a 
fluctuating surface with long-ranged spatial correlation, described either by the KPZ or the EW equations driven by 
a long-ranged noise. Our model systems are semi- autonomous, i.e., we ignore feedback due to the percolating field 
on the environment. The general picture that emerges out of our calculations is that depending upon relative values 
of e = 4 — d and y, a parameter that fixes the spatial scaling of variances of the external forces in the NS, EW or 
KPZ equations, one obtains different universal behavior. However, the details of the ensuing phase diagram in the 
e — y plane depend explicitly on the model used to describe the environment (NS, EW or KPZ). On general ground 
we predict possibilities of four different phases (i) Gaussian, (ii) original DP, (iii) Long-range DP (LDP) and (iv) 
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Long-range (LR). The Gaussian fixed point is generally unstable. The rest are model dependent (i.e., depends upon 
whether the environment is modeled by the NS, EW or KPZ Eq.). For a randomly stirred environment described by 
the NS Eq., we find when e > 4y, the ensuing universal critical behavior of the A APT transition is described by the 
standard DP universality class characterized the DP fixed point with the renormalizcd coupling constants u R — 2e/3, 
w R — 0, corresponding to a dynamic exponent < ze, the dynamic exponent of the environment. Thus one obtains 
weak dynamic scaling, despite the nonlinear coupling between (f> and the environment (here: velocity v). In the 
other regime, i.e., when e < 4y, the DP fixed point gets unstable against perturbations due to the environmental 
fluctuations, and the AAPT is described by a set of critical exponents that depend upon y. The corresponding phases 
are described by either the LDP fixed point u R — 4e/5 — 8y/15 and w R — 8y/15 — 2e/15, or by the LR fixed point 
u R = and w R — 4y/9. In both these regimes = z E = 2 — y/3, thus describing strong dynamic scaling. We 
obtain the relevant scaling exponents in all the phases. Agreement of our results for the scaling exponents and the 
phase diagram with those in Ref. [l8[ shows the general robustness of the perturbation theory, despite our using 
(slightly) different choices for the coupling constants. It also shows how such choices may be exploited to infer the 
large length-scale, long-time limit physics of the system in a simple manner. In contrast, when the environment is a 
fluctuating surface modeled by the EW Eq., the DP phase is linearly stable for e > 3y. We further find that the LR 
phase is not at all stable for any y > 0. Furthermore, the LDP phase does not exist in the range y < e < 7y. Thus, 
the phase diagram in this case does not have any physically meaningful phase in the range y < e < 3y, unlike the NS 
case where the phase diagram is fully spanned by the stable phases of the system. The LDP phase here corresponds 
to strong dynamic scaling and the DP phase weak dynamic scaling. We obtain the associated scaling exponents as 
well. For an environment described by a fluctuating KPZ surface, we again find the existence of four different phases 
similar to the previous cases. We are able to obtain the scaling exponents in each phase and show that strong dynamic 
scaling prevails in the LDP and LR phases, as expected. However, due to the algebraically complicated nature of the 
DRG /3-functions we do not discuss their linear stability here. Validity of our results are limited by the applicability 
of one-loop approximations which suffer from well-known technical difficulties 23[ in systems with long range noises. 
Thus it is important to verify our results in numerical simulations of the models used here. Further, it is unclear 
how multiscaling of the velocity field given by the randomly stirred NS model or the height field given by the KPZ 
equation affect the scaling of the percolating agent <p, or whether <p itself will display multiscaling for its higher order 
structure functions. Since the EW equation is a linear equation, it does not show any multiscaling. However still, field 
(j) may display multiscaling, similar to the passive scalar problem of fluid turbulence, where even if the velocity field 
is Gaussian distributed (albeit with a long ranged spatial correlation), the passive scalar density exhibits multiscaling 
[35j . We look forward to numerical solutions of the continuum model equations in resolving the outstanding theoretical 
issues discussed above. In addition, the effects of EW or KPZ fluctuating surfaces on the DP universality may be 
studied by numerical simulations of lattice-gas type models, which may be constructed by borrowing, e.g., the lattice- 
gas models used in Ref. [36| to study the dynamics of a passive scalar on fluctuating surfaces. Effects of turbulent flow 
on the DP universality may in principle be investigated by adopting the experimental set up of Ref. (ll| and stirring 
the system with long-ranged external forces, or by observing grwoth/decay of a bacteria colony in a turbulent flow. 

Our models and results provide for simple examples of weak dynamic scaling, situations which are not very commonly 
found. Known examples include the studies in Refs. [3Ti - [40l | . Our semi-autonomous models, where effects of the local 
population on the environment is completely neglected, are certainly a simplification of more realistic situations where 
such feedback effects should be present in general. A non-zero feedback is known to affect the scaling properties of 
the NESS in general. For instance, the scaling/multiscaling properties of the magnetic fields in three-dimensional 
Magnetohydrodynamics (3dMHD) are vastly different when the feedback (in the form of Lorentz forces in 3dMHD) 
is present from when it is absent (the passive vector limit) [4l|. Given this, it would be interesting see how feedback 
due to (f> may alter the emerging scaling behavior discussed here. When a feedback is present, the overall system is 
no longer autonomous, but fully coupled. It would be intriguing to see if weak dynamic scaling persists even in the 
fully coupled case. A truly novel feature of our results is that, in all the three models considered here regardless of 
stability issues, the LR and the LDP FPs correspond to the same value of the dynamic exponent z^, but different 
values for the static scaling exponents 77^ , v. This is in contrast with what one finds in equilibrium critical dynamics. 
For instance, model A and model B (in the nomenclature of Ref. 43]) display different dynamic exponent but the 
same static critical exponents for the second order phase transition in the O(N) model. Equality of the static critical 
exponents for different dynamics (i.e., different dynamic exponents) is a requirement of thermal equilibrium. Since our 
work concerns here models that are driven our of equilibrium, such considerations do not arise. Numerical solutions of 
the continuum stochastic equations of motion or simulations of equivalent lattice models should be able to verify our 
results. Beyond immediate theoretical motivation, our work sheds important light for more realistic problems, e.g., 
population dynamics of a bacteria colony or a biofilm resting over a fluctuating surface or a fluctuating biomembrane. 
Further, Our results provide important insight for more biologically motivated problems, e.g., population dynamics 
of a bacteria colony over a fluctuating surface, or in the presence of a macroscopic order, e.g., in a nematic or polar 
active fluid 0. 
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